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Abstract. IWilson fc Bappul l)l957|) have shown the existence of a remarkable correlation between the width of 
the emission in the core of the K line of Ca II and the absolute visual magnitude of late-type stars. 

Here we present a new calibration of the Wilson-Bappu effect based on a sample of 119 nearby stars. We use, 
for the first time, width measurements based on high resolution and high signal to noise ratio CCD spectra and 
absolute visual magnitudes from the Hipparcos database. 

Our primary goal is to investigate the possibility of using the Wilson-Bappu effect to determine accurate 
distances to single stars and groups. 

The result of our calibration fitting of the Wilson-Bappu relationship is My = 33.2 — 18.0 • log Wo, and the 
determination seems free of systematic effects. The root mean square error of the fitting is 0.6 magnitudes. This 
error is mostly accounted for by measurement errors and intrinsic variability of Wo, but in addition a possible 
dependence on the metallicity is found, which becomes clearly noticeable for metallicities below [Fe/H]~-0.4. This 
detection is possible because in our sample [Fe/H] ranges from -1.5 to 0.4. 

The Wilson-Bappu effect can be used confidently for all metallicities not lower than ~-0.4, including the LMC. 
While it does not provide accurate distances to single stars, it is a useful tool to determine accurate distances to 
clusters and aggregates, where a sufficient number of stars can be observed. 

We apply the Wilson-Bappu effect to published data of the open cluster M67; the retrieved distance modulus 
is of 9.65 magnitude, in very good agreement with the best distance estimations for this cluster, based on main 
sequence fitting. 
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1. Introduction 

The K line of Ca II is the deepest and widest absorption line in the spectra of the late type stars, and its core sh ows 
the characteristic double r eversal profile: an emissi on with a central self-absorption (see e.g. IPasauini et aflll988|) . 

Since the discovery by IWilson fc Bappul l|l957t) of the existence of a linear relationship between the logarithm of 
the width of the Call emission (Wo) and the stellar absolute visual magnitude (the so called Wilson-Bappu effect 
), several calibrations of this effect have been attempted. However the Wilson-Bappu relationship (WBR) has only 
seldom been used to determine stellar distances to single stars or aggregates. 

The reliability of past calibrations of the WBR has been limited by the lack of two crucial elements: 

1. Precise and accurate parallaxes for nearby stars; 

2. High quality data and objective method of measuring Wq. 

The first point has been overcome bv IWallerstein et alJ <)l999j) . who produced a new calibration by exploiting the 
Hipparcos database. This work however was hampered by the fact that it was still based on Wq measurements 
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provided bv lWilsonl lll967l ll97fi) , performed on photographic plates. For our calibration we have used Wo measurements 
performed on new, high quality CCD spectra, so that also the second point has been, for the first time, properly 
accounted for. 

A new, reliable WBR determination is especially interesting since new detectors and state of the art spectrographs 
can now produce excellent Call data even for stars in stellar clusters and associations as distant as several kiloparsecs. 
For these clusters and associations, we could therefore apply the WBR to derive their distance. This possible application 
is the main ground of our effort to retrieve a reliable calibration for the WBR. 



2. Data sample and observations 

The full sample for which spectra have been collected (data shown in Tabled consists of 152 stars, but the present 
study is limited to stars with relative parallax errors smaller than 10%. We have also excluded from the original sample 
known multiple systems. After this trimming, the final sample includes 119 stars. 

All the stars but the Sun are included in the Hipparcos catalogue, from which trigonometric parallaxes and visual 
magnitudes have been taken. 

The absolute visual magnitude of the Sun has been taken from lHavesI ()l985|) The sample has been selected in order 
to span a wide range of luminosities (from My — —5 to My ~ 9). 

The observations were obtained between November 1988 and September 1996, at ESO, La Silla, with the Coude 
Echelle Spectrometer, at the focus of the Coude Auxiliary Telescope. The resolution is R=6000 and the S/N ratio 
ranges from ~ 30 to ~ 100 at the bottom of the line (for more details about the first spectra see IPasauinil l|l992lW 

For ~ 30 stars, multiple spectra were taken, and for th e Sun 9 spectra are avai l able. 

Spectral types and metallicities were obtained from the lCavrel de Strobel et aL I <1997l) catalogue. For the stars not 
inclu ded in this cata l ogue, spectral types are taken from the Hipparcos database. The projected rotation velocities are 
from lClebocki et all fcOOOft . 



3. Measurements and calibration 

Wo was measured manually on all of the spectra. We have also developed an IDL macro which performs a multi 
Gaussian fitting to the K-line double-reversal profile, but it has not been used for the final computation because it 
has proven quite fragile for the correct determination of the broad line absorption profile. We are working on improving 
the method, and we may report further progress elsewhere. 

Wq has been computed as the difference in wavelength between the two points taken at the intensity equal to the 
average between those of the Kl minimum and K2 peak on either side of the emission profile (see Figure EJ). 

We found that the definition of Wo which we have adopted correlates better with My than two other widths also 
measured: 

• the difference in wavelength between Kl minima; 

• the width taken at the intensity equal to the average between those of the Kl minima and K2 maxima. 

No correction for instrumental broadening was applied to Wo- The reason is explained in Sect|5.4| 

We notice that our definition of Wo differs slightly from others used in the literature. I Wilson fc Bappul l|l957t) 
define W as the difference in wavelength between the red edge and the violet edge of the emi s sion p rofile. They apply 
to the measured value in the spectra a linear correction: Wo = W — 15 km • s " kept the original 

definition but applied a revised correction: 18 km • s _1 instead of 15. lLut d l|l97Cf) introduced a new definition of Wo, 
defining it as the width at half of the maximum of the emission profile. Lutz's definition is very similar to ours with 
the exception that we found the our definition easier to use in case of difficult spectra and more robust (see below) . 

Wo measurements are subject, of course, to measurement errors. For each star we have computed an accuracy 
qualifier, AWo, in the following way. For the stars with two or more spectra available we derived it as the half of 
the difference between the largest and smallest measured widths. For stars with only one available spectrum, we 
measured Wo with multiple methods, and took the difference between the extrema of the measurements. Then we 
added quadratically to this value (which in some cases was 0) that of the Sun (0.030 A), whose Wo measurements 
variations are supposed to be caused only by the intrinsic variation of the line width and by the limit of the resolution 
power of the spectrograph. We have used A Wo as an estimate of the mean error in Wo. The standard error of log Wo 
(where Wo is in km • s _1 ) for each star, is retrieved applying the propagation of the mean error. 

Independent of the signal to noise ratio, for some stars it is intrinsically more difficult to measure Wo. This is 
because their spectra show asymmetric self absorption, either produced by interstellar lines or by blucshiftcd winds or 
cosmic rays. Inactive, low luminosity stars will typically show shallow reversals, which are more difficult to measure. 

Some of the most doubtful examples and difficult cases are presented in Figure ^ in order to make the reader 
acquainted with the spectra and the possible error sources. In some cases, when strong blending was present or the 
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profile was highly asymmetric, we have measured Wo by doubling the value measured for the "clean" half of the line. 
Anyway, the majority of the spectra we dealt with were as good as the one showed in Figure |21 

We have also computed the standard error on the absolute magnitude of each star. This error has two components: 
the mean error on the apparent visual magnitude (which is, in most of the cases, negligible) and the error given by 
the uncertainty in the parallax, which has to be computed via the propagation of the error. 

The fit of the WB R was performed by means of the IDL routine "fitexy" , which implements the algorithm described 
in lPress et alJ l|l989|) . The algorithm fits a straight line to a set of data points by taking into accounts errors on both 
coordinates. 

We repeated the solution by rejecting stars not passing 3er— , 2.5er- or 2<r-criteria. The results so obtained are 
presented in Table 

It is fundamental to note that, independent of the different sigma clipping criterion used, the solutions found are 
extremely stable, giving the same fit to within la. 

We adopt in the following: My = 33.2 - 18.0 • log W (W is in km • s" 1 ) obtained rejecting HD 63077 and HD 
211998, with a standard deviation: ctwbb, — 0-6 mag. The two rejected stars are the most metal poor of the sample, 
which we will argue in Sect. 15.51 is the most likely cause of their large residuals. 

Hereafter we indicate with My{K) the value of My derived for a single star from its Wq via the WBR. In Figure 
[5] the log Wo vs My diagram is shown, with the error bars representing standard errors in both coordinates. The 
calibration line retrieved is also plotted. 



4. Comparison with other results 

In order to evaluate the external robustness of the relationship, a comparison with other independent investigations 
is fundamental. Many at tempts in the past 40 years have been made to calibrate the WBR based on ground-base d 
trigonometric parallaxes l|Wilsonll959HHodge fc Wallersteirll966tlLutz fe Kelkerll975tlGlebocki fe Stawikowskill978ft . 
Wall erstein et alJ l|l999h (hereafter WMPG) used a very large sample of stars with Hipparcos parallaxes and relative 
measurement errors smalle r than 20%. A s alrea dy stated, their work is based on the width measurements available in 
the literature, mostly from I Wilson! l|l967lll976[l . wich suffer from low spectral acc uracy. 

We note here that such measurements are based on the definition of W given in IWilson fc Bar^nTl (Il957t) . corrected 
for instrumental broadening as in IWilsonl lll959h (see Sect.^J. So the quantity Wo which they adopt is not exactly the 
same as ours, although the two quantities are expected to be strongly correlated. 

WMPG used a linear least squares fitting both not weighted and weighted only in absolute magnitude with As 
WMPG advised, using weighted least squares means giving more weight to the lowest part of the diagram, containing 
the dwarfs that are, on the average, much closer, and therefore with smaller measurement errors on parallaxes. Using 
the weight for both the coordinates, as we did, does not produce to the same effect, because at the same time the 
dwarfs have also smaller Wo, and so larger relative measurement errors, i.e. larger standard errors for log Wo . In Figure 
El we show the comparison of our calibration and the weighted calibration of WMPG. 

We think that this comparison is a strong test of the reliability of the WBR. In fact, with the only exception 
of having used Hipparcos parallaxes, the two samples are independent: less than the 40% of the stars we used in 
the computation are common to WMPG. Different stars means different distributions in magnitudes, metallicities 
and effective temperatures. Furthermore, different spectrographs were used, different measurements (even based on 
different definitions of Wo) were performed as well as a different analysis. Taken at the face values, our calibration and 
WMPG's one are in good agreement for all stars brighter than My ~2, while for dwarfs, the discrepancy is stronger, 
reaching a difference of 1.5 magnitudes at My ~ 9. 

In order to investigate how much of this discrepancy is due to the differences in Wo measurements, we compared, 
for the 64 stars common to the two data sets, the Wo measurements of WMPG and ours. Actually 20 of these stars 
are among the 33 not used in our WBR computation, but for the present comparison this is irrelevant. The two sets 
of measurements show, as expected, a very strong linear correlation: the slope of the Wq out vs Wo wmpg linear fitting 
is very much closed to unity: 1.003, with an intercept of -5.34 km • s _1 (see Figured). 

If we subtract 5.34 km • s _1 to our Wo measurements, we obtain a data set homogeneous to that of WMPG, and 
performing the fitting with the new values gives following WBR: My — 29.7 — 16.3 • log Wo. This result matches very 
well that of WMPG, as it can be seen from Figure We conlcude that the reliability of the WBR is excellent, and 
that the only reason for the discrepancy between the calibrations is the difference in the definition of Wo . 

On the other hand we notice that care is needed in measuring Wo: its definition, and possibly the resolution of the 
spectra used shoud be the same as those of the calibration adopted. 

The fact that the difference between our and WMPG's measurements is about 5 km ■ s , quite similar to the 
projected slit width for R=60000, could suggest that the instrumental profile should indeed be linearly subtracted by 
our Wq measurements to obtain an-instrument free calibration. 

We do not believe that this is the case, because: 
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1. Other studies llLutzll l97m have shown as the case of linear subtraction of instrumental profile is not the best choice. 

2. Our observations of one star (HD 36069) taken at different resolutions to estimate this effect show that the variations 
observed in Wo measurements are not better accounted for by a linear sibtraction of the instrumental profile, as 
shown from the measurements in Table |3J 

For sake of completeness, we remind that the WB R is also valid for the k-line of the Mgll, and that the best 
calibration to date is the one of lCassatella et alJ l|200l|) . which also uses Hipparcos data and IUE spectra. Their result 
is: M v = 34.56 - 16.75 • log W 

5. Is the WBR a good distance indicator? 

The spread around the WBR is still too large to consider it as a reliable distance indicator for single stars. The question 
we are now going to investigate in this section is if the WBR is suitable to determine the distance of clusters of stars. 
A necessary condition which such clusters have to satisfy is, of course, that for a sufficient number of members, a high 
quality spectrum, showing a clear double reversal profile of the K-line, is available. 

The possibility of using the WBR to determine accurate cluster distances is strictly related to the causes of the 
scatter: whether or not it is due to entirely random errors or systematic effects. Among the possible causes of scatter, 
we mention: 

random effects: 

- measurement errors, 

- cyclic variation in the chromospheric activity, 

- variability of some of the stars of the sample; 
systematic effects: 

- reddening, 

- instrumental effects, 

- Lutz-Kelkcr effect (hereafter LKE), 

- hidden parameters, i.e. parameters other than Wo and My on which the WBR could depend. 

White k. Livingstonl ljl98l[) observed the chromospheric emission of the K line of the sun during a whole solar 
cycle. They found a maximum variation of log Wq of about 0.05 during such a period. If we assume that most of 
the stars are affected by a variation of log Wo of the same order of magnitude, the amount of scatter introduced by 
the cyclic variation of the chromospheric activity would represent a relevant fraction of the spread observed in the 
data. Nevertheless, this variability cannot fully explain the observed root mean square error of the WBR fitting. With 
typical uncertainities in Wo due to measurement errors and natural variations of the stellar line width of about 0.036 
A(cf. Table 1), this error, for stars with intermediate widths, say Wq= 0.8 A, accounts for about 0.35 mag of ctwbr 
= 0.6 mag. Therefore, it is necessary to investigate further reasons of uncerta inty in the determinat ion of the WBR. 

Among possible causes of biases we should consider reddening, the LKE l|Lutz fc Kelkeil lll973h . hereafter LKP), 
instrumental effects and the presence of multiple systems. 

5.1. Lutz-Kelker effect 

The LKE is the bias due to the fact that a symmetric error interval [it — ov, it + ov] around the estimated parallax tt, 
does not correspond to a symmetric error interval in distances around — . The inner spherical corona centred in the Sun 
having radii ^ * and i, has a volume smaller than the outer spherical corona. So, assuming a homogeneous space 
density for the stars, we expect that, for a fixed measured parallax, stars having a true distance greater than -, i.e. 
those in the outer corona, will outnumber the stars having a distance smaller than -. There is therefore a systematic 
trend to underestimate distances. The correction which has to be applied to each star, has been calculated in LKP. It 
depends only on the relative error — . 

Our sample has been selected to include only stars with < 0.1). Furthermore, out of the 119 stars, only 7 
have — exceeding 0.075. For these values the LKE is negligible compared with other errors involved: 0.06 mag for 
^ = 0^.075 and 0.11 mag for ^ = 0.1 (See Table 1 in LKP) 

5.2. Reddening 

We have assumed that all stars have zero reddening. This assumption is justified by the fact that the sample stars are 
all within ^200 parsec. 

The most distant star, HD 43455, has a distance of 205 pc, and it is the only one for which we were not able to find 
out a secure upper limit to the reddening. 

HD 78647 has a distance of 176 pc, and a galactic latitude lower than 7.6°, so we can get a rough estimation of its 
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reddening on the Neckel & Klare's maps llNeckel et al.l|l980h. For t his st ar, Ay does not exceed ~ 0.1 mag. 
The remaining stars are within 107 pc. According to Sfeir et all l|l999|) (see their Figure 2) the upper limit of the 
equivalent widt h of the D2 Na I l ine for such a distance is 200 mA. From this quantity we can get the Hydrogen 
column density IjWelsh et al ] |1994 : N(HI)~ 2 • 10 20 , which yields a colour excess: Eg_y ~ 0.03, or an upper limit for 
Ay of about 0.1 mag. 

Fur thermore, 100 of the 119 stars in this sample, are within 75 parsec, so they are in the so called Local Bubble (see 
e.g. lSfeir et al]ll999h . and they are not affected by detectable extinction. 



5.3. Multiple systems 

The presence of unrecognised multiple systems in the sample, gives rise to a systematic underestimation of My(K), 
because we could associate the width of the K line emission of one component to the magnitude of the whole system. 
To avoid this effect, we checked all the objects of our sample on the SIMBAD database and we have excluded all the 
known multiple systems. 



5.4. Instrumental effects 

The measured Wo is likely larger than the intrisic one because of the broadening introduced by the spectrograph. 
The larger the projected slit width is the stronger the intrumental broadening will be. We have shown in Sect. by 
means of data in Tabled that a linear correction for instrumental br oadening; (i.e . subtracting the projected slit width 
from Wo) would not be appropriate. Similar results were found bv iLuta (|1970J), who concluded that a quadratical 
correction should be used. To minimize this effect, our calibration is based on high resolution spectra (the projected 
slit width is about 0.066 Aor 5 km • s _1 ), and appling a quadratical correction even to the smallest Wo value (that of 
HD 42581, 0.30 A) we would obtain: W - W> correctcd = W - - (0.066) 2 = 0.0074A, well below its estimated 

measurement error, i.e. A Wo = 0.02A. For larger values of Wo, Wo — Wo oorrected is even smaller. Hence, the quadratical 
correction is negligible for all stars in our sample. We believe that the quadratical correction is more appropriate than 
the linear one, and it should be applyied when dealing with low resolution spectra, but it is not certain that such a 
small adjustment would represent a real improvement when dealing with data of resolution comparable to that used 
in this work. 



5.5. Other causes 

Many authors searched for ad ditional parameters on which the WBR could depend, finding contradictory results. 
Idebocki fc Stawikowskil l(l978h proposed a corrected WBR with a term for the intensity of the emission that WMPG 
have rejected. 

IP arsons! l)200l|) . analysing the calibration of WMPG, suggested a trend for high luminosity stars that our data 
seem not to confirm: he suggested that O — C (i.e. the difference between the absolute magnitude from Hipparcos 
parallax and the one retrieved by means of the WBR) increases with increasing T e f f for spectral types earlier than ~ 
K3, while the opposite is true for the other stars. He also concludes that this trend gets stronger for brighter stars. 
According to Figurc|Hl while we can draw no conclusions for late type stars, our data seem to suggest a trend opposite 
to that proposed bv lParsonsI (J2001) for spectral types earlier than ~ K3. 

The most obvious hidden parameter to search for is projected rotational velocity. High rotatio nal velocity can 
influe nce Wo in several ways, either because fast rotating stars will tend to be more active (see e.g. ICu^isj^c^.SLal 
2002), or because the width of the line core may be modified by the higher rotational velocity (see e.g. lPasauini et al 
1989). We have 53 stars for which V ■ sini is available, and none are really fast rotators, only for one object V ■ sini 
exceeds 10 km • s _1 . Our conclusion is that, among slow rotators, there is hardly any dependence of the residuals on 
V ■ sin i: we find a correlation coefficient of 0.12 

We have finally searched for a dependence of the O — C on metallicity. Such a dependence can also be expected, 
considering that in stars having lower abundances the core of the line may sample different layers of the atmosphere. 

In particular we have checked whether the WBR is still valid for very metal poor stars. Figure 01 shows two O — C 
vs [Fe/H] diagrams: the one on the left refers to all the stars with available metallicities, in the other diagram only 
stars with [Fe/H]< —0.3 are plotted. A weak but not negligible dependence of the WBR on metallicity does exist, and 
it gets much stronger for metal poor stars. The correlation coefficient is 0.64, and it becomes 0.82 when the 19 most 
metal poor stars are considered, as shown on the right panel of figure 0] 

19 stars are too few to obtain any firm quantitative conclusions. In particular, the O — C vs [Fe/H] relationship, 
to which they would point out (the straight line in right panel of Figure , should be further investigated by means 
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of a richer samp l e. Th e existence of such a relationship for metal poor stars has been independently suggested by 
iDupree fc Sirrithl lll995h . who studied 53 metal poor giants, none of which is in our sample. 

We think that the WBR should be applied very carefully to very metal poor stars (e.g. stars more metal poor 
than [Fe/H]^-0.4) and that further metal poor calibrators should be observed before applying it to very metal poor 
clusters. 



6. Application to M67 

After deriving the WBR, and showing that the scatter is mostly due to random errors, we have the o pportunity to 
test it on a group of stars belonging to a well studied open cluster. M67 Call spectra were published bv lDuoree et alJ 
l)l999|) (See Figure 2, 3 and 4 therein) for 15 stars on the RGB and clump region, and they are suitable for our analysis 
of the WBR. Andrea Dupree kindly provided us with all the spectra in digital form. 

Since M67 has b een extensivel y stud ied, the retrieved distance modulus can be compared with values obtained 
from other authors. ICarraro et alJ (^996) provide a detailed study of M67. They derive, on the basis of the Colour 
M agnitude Diagram, 9 . 55 < (to — M)v < 9.65 mag. 

iMontgomerv et all l)l993(l performed a photometric survey of the central region of M67. They compared their 
photometry with two theoretical isochrones to retrieve distance modulus and age for M67. From V,B — V CMDs, they 
found (m — M)y =9.60 for both i sochro nes (but different ages were found), they have also used a V, V — I CMD, giving 
(m - AfW=9.85. lbinescu et all l|l995|) found 9.7 < (m - M) v < 9.8 mag, obtained by letting E B _ V varying between 
its upper (0.06 mag) and lower (0.0 3 mag:) limits. Their isochrones were constructed using model atmospheres with 
new opacities. In IMontgomerv et alJ l)l993|) other results from the literature are reported, ranging from 9.55 to 9.61. 

In summary, all distance modulus determinations for M67 are in the range 9.55 < (to — M)v < 9.85 mag. 

M67 is a solar metallicity cluster, so we do not need to take care of the metallicity effect which may affect the 
WBR. The Wo measurements were performed in the same way as for the calibration stars, and the results are given 
in Table 0| Out of the 15 stars of the original sample we have selected a subsample of 10, which suitable spectra 
were available, either for quality or clearness of the core reversal. In fact some of the spectra do not show a clear 
unambiguously recognisable double reve rsal feature, s o that the measurement is unreliable. We did not use the stars 
with the following Sanders ID numbers l)Sanderdll977j) : 258, 989, 1074, 1316, 1279. Even among the 10 selected stars 
some show a clearer profile than others, and for four of them the measurements were more uncertain (of the order of 
0.1 A) and they have been flagged with an asterisk in Table 01 

We have to consider that M67 spectra were acquired for other purposes, and in particular they have lower resolution 
and lower S/N ratio than the typical calibration spectra, so we expect a standard error on the single measurement 
higher than the <jwbr derived above. 

In Table^Jthe distance modulus determinations for the single stars retrieved by means of the WBR are given. They 
range from 8.1 to 10.9 mag. The mean value is ~9.7 In spite of the poorer quality of the spectra, all the deviation can 
be explained on the base of the intrinsic spread around the WBR. 

For sake of accuracy we have also taken into account the effect of the difference in resolution between the calibration 
spectra and the M67 observations (5 and 11 km • s _1 respectively). A simple, quadratic correction for the difference 
between the two projected slit widths is applied in the sixth column of Table The correction does not change the 
result in an appreciable way. 

When considering all stars a simple mean gives (M-m)=9.61 mag; which becomes 9.65 when discarding the 4 most 
uncertain measurements. 

We expect that the standard error in our determination of the distance of M 67 would be about: a n ^ R ~ 0.3 if 
we used 6 spectra of quality similar to those used for our calibration (<7' M67 ~ 0.2 if we had 10 spectra of the same 
quality) . 

Trying to push further this application would definitely represent a gross over interpretation of the data, however 
we find it extremely interesting and encouraging that a simple application, using published data, can provide a distance 
modulus in the range between 9.5 and 9.8, in excellent agreement with completely independent measurements, such 
as those obtained with main sequence model fitting. 

7. Conclusions 

We have shown that the coupling of CCD high resolution, high S/N ratio data with the use of the Hipparcos parallaxes 
allows a good determination of the WBR. The root mean square error found around this relationship (0.6 magnitude) 
is not good enough to determine accurate distances to single stars, but it can be used to infer accurate distances 
of clusters or groups, provided that they are not too metal poor. This is possible because the uncertainties in the 
relationship are mostly due to random errors (measurements, cycles) and not from systematic effects. This implies 
that once one has observed a sufficient number of stars, n, the distance modulus standard error can be reduced to about 
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O.Gmag/v 7 ^- Its extension to metal poor objects (e.g. stars with Fc/H <-0.4) would require extra care to fully evaluate 
the impact of low metallicity on the relationship. When using our WBR in photometric parallax determinations, the 
resolution used should be comparable (within a factor <~ 3) to that of the calibration (R=60000), to avoid large 
corrections, and care has to be exercised in measuring Wo, following the proper calibration definition. 
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the referee, Elena Schilbach, for very valuable comments and suggestions, which improved considerably the quality of this paper. 
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Tabic 1: Our data sample. The full sample consists of 152 stars, 
33 of which have not been used in the present study, because of 
their high measurement error in the parallax (— > 0.1) or because 
they are multiple systems. We have flagged their HD identificators 
with an asterisk. Column 1: HD ID number of the star. Column 
2: Wq in A. Column 3: accuracy qualifier for Wo, again in A. Its 
mean value is 0.036. Column 4: absolute visual magnitude using 
Hippa rcos parallaxes. For the Sun we used the value given bv lHavea 
(1985). Column 5: spectral type from Cayrel De Strobel Catalogue. 
When the spectral type is not available from this catalogue the data 
is taken from the Hipparcos database, and is flagged with a dagger. 
Column 6: mean metallicity, when available from Cayrel De Strobel 
Cata logue. Column 7: proj ected rotation velocity, when available 
from IClebocki et all ll2000l) . 



Star 


W [A] 


AWo [A] 


M v 


Sp.Type 


[Fe/H] 


V ■ sin i 


SUN 


0.49 


0.030 


4.82 


G2V 


0.00 


1.6 


HD 203244 


0.47 


0.030 


5.42 


G5V 


-0.21 


- 


HD 17051 


0.62 


0.030 


4.22 


GOV 


-0.04 


5.7 


HD 1273* 


0.5 


0.030 


5.03 


G2V f 


-0.61 


- 


HD 20407 


0.41 


0.042 


4.82 


G1V 


-0.55 


- 


HD 20766 


0.50 


0.030 


5.11 


G2.5V 


-0.25 


- 


HD 20630 


0.53 


0.030 


5.03 


G5Vvar 


0.11 


4.6 


HD 20807 


0.44 


0.036 


4.83 


G1V 


-0.21 


- 


HD 20794 


0.43 


0.030 


5.35 


G8V 


-0.38 


- 


HD 26491 


0.50 


0.058 


4.54 


G3V 


-0.23 


- 


HD 1581 


0.49 


0.15 


4.56 


F9V 


-0.20 


3 


HD 30495 


0.55 


0.030 


4.87 


G3V 


0.11 


3 


HD 32778 


0.42 


0.030 


5.28 


GOV 


-0.61 


- 


HD 34721 


0.60 


0.032 


3.98 


GOV 


-0.25 


- 


HD 36435 


0.48 


0.030 


5.53 


G6-G8V 


-0.02 


4.5 


HD 39587 


0.56 


0.030 


4.70 


GOV 


0.08 


9.3 


HD 43834 


0.51 


0.036 


5.05 


G6V 


0.01 


1.8 


HD 48938 


0.53 


0.042 


4.31 


G2V 


-0.47 


- 


HD 3443 


0.47 


0.036 


4.61 


K1V 


-0.16 


2.7 


HD 53705 


0.60 


0.036 


4.51 


G3V 


-0.30 


- 


HD 3795 


0.46 


0.030 


3.86 


G3-G5V 


-0.73 


- 


HD 63077 


0.41 


0.032 


4.45 


GOV 


-0.90 


- 


TTr\ onfio 

HD 3823 


0.60 


0.042 


3.86 


G1V 


-0.35 


3 


HD 64096* 


0.50 


0.030 


4.05 


G2V f 






HD 65907 


0.52 


0.050 


4.54 


GOV 


-0.36 




HD 67458 


0.47 


0.036 


4.76 


G4IV-V 


-0.24 




HD 74772 


0.89 


0.032 


-0.17 


G5III 


-0.03 


5.8 


HD 202457 


0.61 


0.036 


4.13 


G5V 


-0.14 




HD 202560 


0.35 


0.030 


8.71 


M1-M2V f 






HD 202628 


0.55 


0.030 


4.87 


G2V 


-0.14 




HD 202940* 


0.6 


0.030 


5.20 


G5V 


-0.38 


1.2 


HD 211415 


0.48 


0.032 


4.69 


G3V 


-0.36 


1.7 


HD 211998 


0.47 


0.036 


2.98 


A3V 


-1.50 




HD 212330 


0.51 


0.067 


3.75 


G3IV 


0.14 


1.8 


HD 212698 


0.54 


0.030 


4.04 


G3V 


0.08 


9.7 


HD 14412 


0.40 


0.042 


5.81 


G5V 


-0.53 




HD 14802 


0.62 


0.032 


3.48 


G2V 


0.10 


3 


HD 104304 


0.55 


0.036 


4.99 


G9IV 


0.17 


1.7 


HD 114613 


0.56 


0.030 


3.29 


G3V f 




2.7 


HD 11695 


1.02 


0.030 


-0.57 


M4III f 






HD 194640 


0.46 


0.032 


5.17 


G6-G8V f 
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Star 


Wo 


[A] 


AW [A] 


M v 


Sp.Type 


[Fe/H] 


V ■ sin i 


HD 20610* 





,86 


0.032 


0.39 


KOIII 


-0.07 


- 


HD 209100 


0. 


.39 


0.030 


6.89 


K4.5V 


0.14 


0.7 


HD 211038 


0. 


,56 


0.032 


3.64 


K0-K1V 


- 


- 


HD 219215 


1. 


.03 


0.030 


0.05 


M2III f 


- 


- 


HD 29503* 


0. 


.80 


0.030 


1.23 


KOIII 


-0.11 


- 


HD 35162 


0. 


.86 


0.042 


0.28 


G8-K0II-III 


-0.31 


- 


HD 36079* 


0. 


.92 


0.030 


-0.63 


G5II 


-0.20 


5.1 


HD 4128 


0. 


.94 


0.032 


-0.30 


KOIII 


-0.01 


3.3 


HD 43455 


1. 


.05 


0.032 


-1.55 


M2.5III f 


- 


- 


HD 4398 


0. 


.82 


0.032 


0.44 


G8-K0III f 


- 


- 


HD 102212 


1. 


.01 


0.067 


-0.87 


M0III f 


- 


- 


HD 111028 


0. 


.72 


0.036 


2.40 


K1III-IV 


-0.40 


1.5 


HD 112300 


1. 


.04 


0.030 


-0.57 


M3III 


-0.09 


- 


HD 113226 


0. 


.94 


0.030 


0.37 


G8IIIvar 


0.04 


2.8 


HD 114038 


0. 


.89 


0.030 


0.29 


Klin 


-0.04 


- 


HD 115202 


0. 


.66 


0.030 


2.26 


K1III f 


- 


- 


HD 115659 


0. 


.92 


0.030 


-0.04 


G8III 


-0.03 


4.2 


HD 117818 


0. 


.81 


0.032 


0.67 


KOIII 


-0.40 


- 


HD 119149* 


1. 


.15 


0.030 


-0.70 


M2III f 


- 


- 


HD 120477 


0. 


.90 


0.032 


-0.33 


K5.5IIIvar 


-0.23 


2.2 


HD 121299 


0. 


.83 


0.032 


0.70 


K2III 


-0.03 


- 


HD 123123 


0. 


.83 


0.030 


0.79 


K2III f 


-0.05 


- 


HD 124294 


0. 


.95 


0.030 


0.00 


K2.5IIIb 


-0.45 


- 


HD 125454 


0. 


.82 


0.030 


0.52 


G8III 


-0.22 


- 


HD 126868 


0. 


.77 


0.030 


1.72 


G2III f 


- 


14.3 


HD 129312* 


1. 


.04 


0.030 


-1.37 


G7IIvar 


-0.30 


6.5 


HD 130952 


0. 


.88 


0.030 


0.82 


G8III 


-0.29 


- 


HD 133165 


0. 


.83 


0.032 


0.64 


K0.5IIIb 


-0.22 


- 


HD 136514 


0. 


.85 


0.030 


0.98 


K3IIIvar 


-0.14 


0.6 


HD 138716 


0. 


.68 


0.032 


2.30 


K1IV 


-0.13 


2.5 


HD 140573 


0. 


.88 


0.030 


0.87 


K2IIIb 


0.14 


1.4 


HD 141680 


0. 


.82 


0.032 


0.68 


G8III 


-0.28 


1.1 


HD 145001* 


0. 


.94 


0.050 


-0.37 


G8III 


-0.26 


9.9 


HD 145206* 


1. 


.00 


0.032 


-0.51 


K4III 


0.04 


3.2 


HD 146051 


1. 


.12 


0.030 


-0.85 


M0.5III 


0.32 


- 


HD 146791 


0. 


.83 


0.030 


0.64 


G9.5IIIb 


-0.13 


- 


HD 148349* 


0. 


.90 


0.030 


-0.71 


M2 f 


- 


- 


HD 148513* 


1. 


.10 


0.036 


-0.15 


K4III 


-0.14 


0.6 


HD 150416* 


1. 


.04 


0.030 


-0.48 


G8II III 


+0.04 


- 


HD 151217 


1. 


.00 


0.032 


0.00 


K5IIIvar 


-0.11 


2.3 


HD 152334 


0. 


.98 


0.030 


0.30 


K4III f 


- 


- 


HD 152601 


0. 


.82 


0.030 


0.82 


K2III 


0.00 


- 


HD 161096 


0. 


.91 


0.030 


0.76 


K2III 


0.05 


2.7 


HD 164349* 


1. 


.17 


0.030 


-1.84 


K0.5IIb 


-0.32 


- 


HD 165760 


0. 


.91 


0.095 


0.33 


G8III 


-0.15 


2.2 


HD 168723 


0. 


.67 


0.030 


1.84 


K0III-IV 


-0.10 


2.6 


HD 169156* 


0. 


.77 


0.030 


0.82 


G9IIIb 


-0.17 


- 


HD 169767 


0. 


.74 


0.030 


1.15 


G8-K0III f 


- 


- 


HD 170493 


0. 


.44 


0.030 


6.67 


K3V f 




3.5 


HD 171443 


0. 


.98 


0.030 


0.21 


K3III 


0.09 


1.8 


HD 171967* 


1. 


.10 


0.032 


-1.57 


M2III f 






HD 173009* 


1. 


.11 


0.032 


-1.14 


G8IIb 


0.05 


6.0 


HD 173764* 


1. 


.80 


0.036 


-2.40 


G4IIa 


-0.15 


6.5 


HD 175775* 


1. 


.07 


0.030 


-1.76 


G8-K0II-III 


-0.19 
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Star 


Wo [A] 


AW [A] 


M v 


Sp.Type 


[Fe/H] 


V ■ sin i 


HD 176678 


0.84 


0.032 


0.73 


KHIIvar 


-0.19 


- 


HD 177565 


0.51 


0.030 


4.98 


G5IV 


0.03 


- 


HD 17970 


0.44 


0.050 


6.01 


K1V f 


- 


- 


HD 181391* 


0.77 


0.030 


1.61 


G8III 


-0.21 


2.8 


HD 182572 


0.60 


0.032 


4.27 


G8IV 


0.38 


2.3 


HD 183630* 


1.02 


0.050 


-0.90 


MHIvar f 


- 


- 


HD 184406 


0.79 


0.030 


1.80 


K3IIIb 


0.05 


1.3 


HD 186791* 


1.32 


0.030 


-3.02 


K3II 


0.00 


3.5 


HD 188310 


0.91 


0.032 


0.73 


G9IIIb 


-0.32 


2 


HD 189319 


1.13 


0.030 


-1.11 


K5III f 


- 


- 


HD 190248 


0.54 


0.030 


4.62 


G7IV 


-0.26 


- 


HD 190406 


0.49 


0.030 


4.56 


G1V f 


- 


5 


HD 191408 


0.38 


0.030 


6.41 


K2V f 


- 


0.1 


HD 194013 


0.75 


0.032 


0.91 


G8III-IV 


-0.03 


1 


HD 195135 


0.88 


0.036 


1.07 


K2III 


0.03 


- 


HD 196574* 


0.86 


0.030 


-1.04 


G8III 


-0.13 


3.7 


HD 196758 


0.89 


0.036 


0.77 


K1III 


-0.12 


1.8 


HD 196761 


0.42 


0.030 


5.53 


G8-K0V f 


- 


- 


HD 198026* 


1.03 


0.030 


-1.24 


M3IIIvar f 


- 


- 


HD 201381 


0.73 


0.032 


1.00 


G8III 


-0.15 


2.8 


HD 203504 


0.86 


0.030 


0.71 


K1III 


-0.14 


1.2 


HD 205390 


0.41 


0.030 


6.30 


K2V f 


- 


- 


HD 206067 


0.87 


0.032 


0.76 


KOIII 


-0.17 


1 


HD 206453 


0.82 


0.030 


-0.03 


G8III 


-0.20 


- 


HD 206778* 


1.78 


0.030 


-4.19 


K2Ibvar 


-0.05 


6.5 


HD 209747 


1.06 


0.030 


0.32 


K4III 


0.00 


2.3 


HD 209750* 


2.14 


0.030 


-3.88 


G2Ib 


0.18 


6.7 


HD 211931* 


0.86 


0.030 


1.03 


A1V f 


- 


- 


HD 212943 


0.78 


0.030 


1.33 


KOIII 


-0.33 


0.6 


HD 213042 


0.44 


0.030 


6.71 


K4V f 


- 


- 


HD 2151 


0.62 


0.030 


3.45 


G2IVvar 


-0.18 


3 


HD 216032* 


1.03 


0.030 


-1.28 


K5II f 


- 


- 


HD 217357 


0.33 


0.030 


8.33 


K5-M0V f 


- 


- 


HD 21749 


0.38 


0.030 


7.01 


K5V f 


- 


- 


HD 217580 


0.45 


0.036 


6.34 


K4V j 


- 


3.6 


HD 218329 


1.12 


0.030 


-0.43 


M2III f 


- 


- 


HD 220339 


0.42 


0.030 


6.35 


K2V f 


- 


5.5 


HD 220954 


0.86 


0.030 


0.83 


K1III 


-0.10 


0.6 


HD 27274 


0.41 


0.030 


7.06 


K5V f 


- 


- 


HD 32450 


0.33 


0.030 


8.66 


MOV f 


- 


- 


HD 42581 


0.30 


0.030 


9.34 


M1-M2V f 


- 


3 


HD 4747 


0.46 


0.030 


5.78 


G8-K0V f 


- 


- 


HD 56855* 


1.76 


0.030 


-4.91 


K3Ib f 


- 


- 


HD 59717* 


1.05 


0.030 


-0.50 


K5III f 


- 


- 


HD 68290* 


0.84 


0.030 


0.95 


KOIII 


-0.03 


- 


HD 73840* 


1.01 


0.030 


-0.56 


K3III 


-0.21 


- 


HD 74918* 


0.73 


0.030 


0.11 


G8III 


-0.20 


- 


HD 75691 


0.96 


0.030 


-0.01 


K3III 


-0.11 


- 


HD 78647 


1.65 


0.030 


-3.99 


K4Ib-IIvar 


0.23 


8.9 


HD 81101 


0.82 


0.030 


0.62 


G6III f 






HD 82668 


1.21 


0.361 


-1.15 


K5III j 






HD 85444 


0.93 


0.030 


-0.50 


G6-G8III 


-0.14 


2.9 


HD 90432 


1.04 


0.030 


-1.15 


K4III 


-0.12 




HD 93813 


0.96 


0.030 


-0.03 


K0-K1II 


-0.32 
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Star 


W [A] 


AW [A] 


M v 


Sp.Type 


[Fe/H] V-sini 


HD 95272 


0.90 


0.030 


0.44 


KlUI 


-0.15 


HD 9540 


0.50 


0.030 


5.52 


KOV f 




HD 98430 


0.89 


0.030 


-0.31 


KOIII 


-0.40 1.8 
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criterion 


NUMBER OF STARS 


NUMBER OF 


FINAL RESULT 




REJECTED 


ITERATIONS 


a 


b 


(7 WBR 


3(7 


2 


2 


33.2 


-18.0 


0.60 


2.5a 


5 


3 


33.5 


-18.3 


0.56 


2(7 


16 


9 


33.6 


-18.3 


0.49 








<7 a = 


0.5 & 


a b = 0.3. 



Table 2. The resulting calibrations after three different sigma clipping criteria. <jwbr is the the standard deviation 
of a single measurement. a a and are the standard deviations of the retrieved parameters, respectively, a and b of 
My = a + b ■ log WV The stability of the solution shows the reliability of the WBR as a distance indicator. 



Resolution Wo [A] W 0lmcorr [K] W %uadcorr [k] 



R= 110000 


0.923 


0.887 


0.922 


R=80000 


0.915 


0.866 


0.914 


R=60000 


0.925 


0.859 


0.923 


R=60000 


0.918 


0.852 


0.916 


R=40000 


0.929 


0.831 


0.924 


R=30000 


0.953 


0.822 


0.944 


Variance 


1.8 ■ 10 -4 


5.6 • 10" 4 


1.2 ■ 10" 



Table 3. Six Wg measurements of HD 36079, made on spectra with different resolutions, both corrected and not 
corrected for instrumental broadening subtracting linearly the projected slit width. The variance of data in Column 3 
data is about 3 times that of Column 1 and, most important, the corrected Wq decreases with increasing correction. 
These measurements show that the linear subtraction of the instrumental profile is not appropriate for our data. A 
quadratic correction seems to be more justified for our data (see the last Column). 



Column 1 


Column 2 


Column 3 


Column 4 


Column 5 


Column 6 


Sanders ID 


Wo (A) 


my 


Mv{K) 


(m-M)v 


(m-M)v 








no corr 


no corr 


corr 


S1010 


0.851 


10.48 


0.657 


9.823 


9.742 


S1016* 


0.700 


10.30 


2.184 


8.116 


7.996 


S1074* 


0.784 


10.59 


1.298 


9.292 


9.196 


S1135 


0.963 


9.37 


-0.310 


9.680 


9.617 


S1221 


0.854 


10.76 


0.629 


10.131 


10.050 


S1250* 


0.997 


9.69 


-0.581 


10.271 


10.212 


S1479 


0.868 


10.55 


0.502 


10.048 


9.970 


S1553 


0.970 


8.74 


-0.366 


9.106 


9.044 


S488 


1.010 


8.86 


-0.682 


9.542 


9.485 


S978* 


1.080 


9.72 


-1.206 


10.926 


10.876 




Mean value using afl stars: 


9.693 


9.619 


Mean value using only unffagged stars: 


9.722 


9.651 



Table 4. Data about the sample of the 10 stars in M67. Column 1: Sanders ID number l|Sanderall977j) of the star. 
The stars with doubtful measurements are flagged with an ast erisk. Column 2: Wilson-Bappu width in A. Column 3: 
apparent visual magnitude, from Table 1 in lDupree et alJ l)l999|) (see references therein). Column 4: absolute magnitude 
inferred from the WBR. Column 5: retrieved distance modulus. Column 6: retrieved distance modulus using corrected 
widths: W 0corrected = ^/W$ - ^SWjUsoooo + P5W fl= 6 oooo ( PSW is the projected slit width). 



14 



G. Pace et al.: Wilson-Bappu effect 





Fig. 1. Doubtful examples of Call K line profiles. The spectra are affected by cosmic rays, may be blended with 
interstellar absorption, which strongly influence the emission profile observed in HD 82668, or, as in the case of HD 
i02212, show a blucshiftcd wind. In other cases, such as in HD 1581, the emission is very weak. 
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A 

Fig. 2. Spectrum of HD 4128. Most of the spectra of our sample have a comparable quality. 
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Fig. 3. The spectral type diagram vs — C. The spectral types are indexed as in lParsonsl ( 200lh : is for GO stars, 1 
for Gl and so on. Only luminous stars are plotted (no IV and V luminosity classes). 
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Fig. 4. The O — C vs [Fe/H] diagrams both for all stars of the sample with available metallicities (on the left) and for 
metal poor stars only. On the latter is also shown the retrieved regression line, which has a slope as high as 2.61. 




Fig. 5. Our calibration of the Wilson-Bappu Effect: My = 33.2 — f8.0 • log Wo- This calibration is the 3cr 
criterion one in Table [21 (HD 63077 and HD 211998 are not used). The error bars represent standard errors in both 
the coordinates. 
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Fig. 6. Comparison between the following calibrations: the present paper's one [My = 33.2 — 18.0 • log Wo), 
the one we obtained after subtracting 5.34 km • s _1 to all Wo measurements and that of WMPG (the weighted one: 
My — 28.83 — 15.82 • log Wo). The points refer to our data, including stars not used in the calibration because 
of the uncertainty in the parallax > 0.1) or because they are binary or multiple systems. 
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Fig. 7. Comparison between our measurements of Wq and those used by WMPG. 
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